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ANALYSIS OF HELICOPTER ROTOR BLADE 


TORSIONAL OSCILLATIONS DUE TO STALL 

By Peter Crimi 
Avco Systems Division 


SUMMARY 


An analysis of stall -induced helicopter rotor blade torsional oscillations 
was carried out, the primary objectives being to predict the onset and severity 
of the oscillations and their relationship to aircraft and blade parameters. 

Blade flapping, flapwise bending and torsional degrees of freedom were taken 
into account, with radial variation in aerodynamic loading determined from a 
previously developed dynamic stall model, using a strip approximation. 

Results of analyses were compared with data from flight tests of two 
different helicopters. The severity of torsional oscillations and their relation 
to aircraft forward speed agreed fairly well, while the radial and azimuthal 
extent of stall and the azimuthal variation of section loading were in good 
agreement. 

Analyses were carried out while parametrically varying blade vibra- 
tional characteristics. It was found that the amplitudes of the higher harmonics 
of torsional oscillations can be significantly reduced by either reducing the 
torsional natural frequency or introducing viscous damping in the torsional 
degree of freedom. 

A preliminary investigation was also carried out to determine the 
feasibility and practicality of alleviating the stall problem by means of 
boundary layer control. The results indicate that boundary layer control 
would be effective in reducing the higher harmonics of torsional oscillations 
due to stall and that its implementation would not require excessive power or 
suction rates. 



INTRODUCTION 


A problem currently of particular concern in helicopter design and 
operation is the occurrence of large -amplitude oscillatory blade torsional 
deflections and control -linkage loads caused by stall at high forward speed. 

The oscillations effectively limit the speed of the aircraft because the control 
linkages are not able to withstand sustained oscillatory loading of that 
severity (Ref. 1). 

Numerous studies of various aspects of the problem have been conducted 
in the past several years. Many of these have been concerned with the 
process of unsteady airfoil stall and its relationship to excitation of the tor- 
sional oscillations. Experimental studies (Refs. 2-5, for example) have 
shown that the dominant feature of dynamic stall is the development and 
growth of a region of vortical flow, starting from the leading edge, which 
induces lift considerably in excess of what is achieved statically and, of 
particular importance to the rotor problem, causes a rapid rearward shift in 
the center of pressure, resulting in severe transient nose -down moments. 
Analyses of unsteady stall of two-dimensional airfoils with viscous effects 
taken into account empirically have been performed (Refs. 6 and 7). More 
recently, an analytical model of unsteady stall was developed which accounts 
for interaction of the viscous and inviscid flows (Refs. 8 and 9). This theory, 
which shows quite good qualitative and fair quantitative agreement with 
measurements of airfoil loading, was employed in the study reported here, 
and is described briefly in the next section. 

Considerable research has also been directed specifically to the rotor 
problem. Results of wind tunnel and flight tests of rotors undergoing stall- 
induced torsional oscillations are reported in Refs. 10 through 13. Analytical 
studies utilizing empirical aerodynamic representations are presented in 
Refs. 14 and 15. The analyses are for the most part in reasonably good 
agreement with the tests, but the use of empirical methods restricts their 
applicability. The study reported in Ref. 16 employed the analytic dynamic 
stall model of Ref. 8, but a two-dimensional elastomechanical system was 
used to approximate that of a rotor blade, the objectives being primarily to 
qualitatively evaluate and determine the interrelationships of different rotor 
aeroelastic instabilities. 

This study was specifically directed to predicting rotor blade torsional 
oscillations due to stall and their relationship to flight conditions and rotor 
parameters, using a more realistic aeroelastic model which includes an 
analytic representation of unsteady stall and accounts for radial variations 
in aerodynamic loading. Responses of the blade in flapping, flapwise bending 
and torsion have been analyzed by forward integration in time with the 



aerodynamic loading evaluated at either five or six radial stations, using the 
strip approximation and assuming uniform inflow. Results of analyses of 
specific rotor configurations have been compared with flight test data, and 
the effect of parametric variations in blade vibrational characteristics 
examined. In addition, a preliminary assessment has been made of the 
possible, use of boundary layer control in the suppression of torsional oscilla- 
tions. 


SYMBOLS 


blade semichord, m 

blade section lift coefficient, ratio of lift per unit span to 
P U 2 b 

blade section moment coefficient, ratio of moment per unit span 
about quarter -chord, positive nose up, to 2 p b 

pressure coefficient, ratio of pressure, less free-stream 
pressure, to p U^/2 

suction coefficient, Cq - Q/(Uc) 

2 4 

rotor thrust coefficient, ratio of thrust to tt p Q R 
blade chord, m 

first coupled torsion-control deflection mode shape 

first flapwise bending mode shape 

average blade torsional stiffness, N-m 

section plunging rate, m/s 

ratio of section moment of inertia to Mg R 

control stiffness, N-m/rad 

3 

mass ratio, k m - p R / 

pushrod moment arm with respect to pitch axis, m 
length of controlled blade segment, m 
blade mass, kg 

blade mass per unit span, kg/m 

dimensionless mass per unit span, m = R Mg/Mg 
number of blades 

total power required by boundary layer control system, hp 
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n 

suction rate per unit span, m^/s 
total suction rate per blade, m /s 

magnitude of fluid velocity at the edge of the boundary layer at 
the separation point, m/s 

rotor radius, m 

Reynolds number, Re = Uc/z/ . 

blade-fixed radial coordinate, m 

inner radius of controlled segment, m 

outer radius of controlled segment, m 

radius at center of controlled segment, m 

ratio of blade root cut-out radius to R 

dimensionless radial coordinate, s = r/R 

magnitude of free -stream velocity relative to blade section, m/ s 

helicopter forward speed, m/s 

helicopter gross weight, N 

wake -induced inflow, m/s 

shaft -fixed rotor coordinates, m 

chordal coordinate with origin at the leading edge, positive aft, m 
ratio of blade section mass center offset, positive aft, to R 
location of downstream end of trapped-air region 
location of separation point 
angle of attack, deg 

shaft angle with respect to vertical, positive aft, rad 
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flapping angle, rad 

possible increase in C across leading-edge bubble 
required increase in C across leading-edge bubble 
ratio of flapping hinge offset to R 

viscous damping coefficient in torsion, fraction of critical 
angle between blade section chord line and rotor plane, rad 
blade root pitch angle, rad 

blade tip angular deflection due to first coupled torsion-control 
deflection mode, rad 

advance ratio, ~ ^f/{ ^ R) 

kinematic viscosity of air, m 4 /s 

1/2 

stall parameter, “ (1 “ Ap/ A r ) 

•3 

air density, kg/m 

blade tip deflection due to first flapwise bending mode, m 
blade azimuth angle, deg or rad 
rotor angular speed, rad/sec 

ratio of nonrotating bending vibration frequency to ft 

ratio of rotating bending vibration frequency to ft 

ratio of nonrotating torsion-control deflection vibration 
frequency to ft 

ratio of rotating torsion-control deflection vibration frequency 
to ft 
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FORMULATIONS 


System Analyzed 

The motions of the rotor blade are referred to shaft -fixed coordinates 
(x, y, z), with the x-axis directed downwind, as indicated in Fig. 1. Blade 
deflection due to flapping and flapwise bending at radial position r is 
/3r + f (p(x) 0, where /3 is flapping angle, 0 is bending deflection at the 
tip (r = R) and f^ is the mode shape of the first flapwise bending mode 
(see Fig. 1). 

The flow relative to a blade section at radial station r, with the blade 
at azimuth angle iff , is as shown in Fig. 2. The plunging rate h indicated 
in the figure results from blade flapping and flapwise bending, and is given by 


h 

OR 



+ fi ( ft + f '<£ d> ) cos \Jj 


The blade section angle 0, measured with respect to the x-y plane, is the 
sum of the built-in twist f Q (r), blade root pitch setting © Q (if/) and tor- 
sional deflection 0,, being the torsional deflection at the tip and Iq 

the mode shape of the first coupled torsion-control deflection mode. The 
section pitch rate is given by 


e dd o d0 i 

n • If * '« <r) 77 * P * '* + 


The last two terms in the above expression derive from the component of 
rotor angular velocity parallel to the blade elastic axis. 

Aerodynamic Loading 

The radial distribution of blade section lift and moment were obtained 
by applying the theory of Refs. 8 and 9 stripwise at up to six radial locations. 
It was necessary to restrict the analysis by only allowing for the occurrence 
of leading-edge stall, even though the theory can take account of trailing -edge 
stall as well, in order to keep computer storage and time requirements within 
practical bounds. Since, as discussed further in the next section, the 
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computed loading does agree quite well with flight test data, this simplification 
apparently was a reasonable one. 

The flow elements represented in the analysis are sketched in Fig. 3. 
When the airfoil is not stalled (Fig. 3a) these are: (1) the laminar boundary 

layer from the stagnation point to separation near the leading-edge, (2) the 
small leading-edge separation bubble, and (3) a potential flow, including a 
vortex wake generated by the variation with time of the circulation about the 
airfoil. When the airfoil is stalled, as indicated in Fig. 3b, the flow elements 
are: (1) the laminar boundary layer, (2) a trapped-air region extending from 

the separation point to the pressure recovery point, and (3) a potential flow 
external to the airfoil and trapped-air region, again including a vortex wake. 
The analytic representations of these elements are described briefly below. 
Details are given in Refs. 8 and 9. 

Potential Flow. — Given the airfoil section characteristics and motions, 
together with the distribution of pressure in the trapped-air region if the air- 
foil is stalled, the flow and pressure over the airfoil must be determined to 
compute the integrated load and analyze the boundary layer. The problem was 
formulated by imposing linearized boundary conditions of flow tangency and 
pressure, using a perturbation velocity potential derived from source and 
vortex distributions. The resulting coupled set of singular integral equations 
is solved by casting the singularity distributions in series form and solving 
for the unknown coefficients by imposing boundary conditions at prescribed 
points. A consistent second-order correction to the linearized solution is 
then calculated and added to the latter solution. 

Boundary Payer. — The boundary layer representation of Ref. 8 employs 
the method of finite differences for unsteady flows with variable step size in 
both streamwise and normal directions, with the error in each finite -differ ence 
approximation the order of the square of the step size. It was determined 
from preliminary calculations performed for the study reported in Ref. 16 
that, at least for leading-edge stall, results are virtually unaffected by assum- 
ing quasi-steady flow in the boundary layer. That assumption was therefore 
employed for all computations, to take advantage of the resulting substantial 
savings in computer storage requirements and computing time. 

Trapped-Air Region. — The function of the model of the trapped-air 
region is to define the streamwise distribution of pressure in that region, 
given the locations of the separation and recovery points and the pressure 
at the recovery point. The region is assumed to consist of a laminar constant- 
pressure free shear layer from separation to transition, a turbulent constant- 
pressure mixing region, and a turbulent pressure -recovery region. The 
laminar shear layer is analyzed by the method of Ref. 17, assuming quasi- 
steady flow. The turbulent mixing and pres sure -recovery regions are 
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(1) LAMINAR BOUNDARY LAYER 



(a) ATTACHED FLOW 


<1) LAMINAR BOUNDARY LAYER 

(2) TRAPPED-AIR REGION 

(3) AIRFOIL AND VORTEX WAKE 



(b) LEADING-EDGE STALL 


Figure 3 FLOW ELEMENTS 
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analyzed using the steady -flow momentum integral and first moment equations. 
Profile parameters in these regions are assumed to be universal functions of 
a dimensionless streamwise coordinate, with those functions derived from an 
exact viscous -inviscid interaction calculation. Matching of approximate 
solutions for the mixing and pressure -recovery regions at their interface 
completes the analysis. 

Leading-Edge Bubble. — The leading -edge bubble on an unstalled airfoil 
is analyzed using the same basic relations employed for the trapped-air 
region. Given the boundary -layer parameters at separation, the length of the 
bubble and the amount of pressure rise possible, for that length, in the pres- 
sure recovery region, are computed. That pressure rise is compared with 
the rise in pressure in the potential flow over the length of the bubble. If the 
latter is greater than the former, the bubble is assumed to have burst, and 
the stall process is initiated. 

Loading Calculation Procedure . — Calculation of the loading at a given 
blade section proceeds by forward integration in time, using the blade motions 
derived by integrating the equations of motion of the elastomechanical system. 
If, at a given instant, the blade section is not stalled, the potential flow is 
computed, and the boundary layer and leading -edge bubble are analyzed to 
check for bubble bursting. If the section is stalled, the pressure distribution 
in the trapped-air region is computed, the potential flow evaluated, and the 
boundary layer is analyzed to locate the separation point. The last two steps 
are repeated iteratively until assumed and computed separation points agree. 
Rate of growth of the trapped-air region is determined from an estimate of 
the rate of fluid entrainment derived from the potential -flow solution. Unstall 
is determined by first postulating its occurrence and analyzing the leading- 
edge bubble which would then form to ascertain whether that event did in fact 
occur. 


Elastomechanical Representation 

As noted previously, a blade with flapping, flapwise bending and torsional 
degrees of freedom was analyzed. Using blade azimuth angle ft as indepen- 
dent variable, the equations of motion for the blade can be written in the fol- 
lowing dimensionless form (Ref. 1 6 ) : 



+ 





i = I, 2, 3; 
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where = P , = 0/R and Z 3 = The ^Iq's are nonrotating, un- 

coupled natural frequencies ( <0 l o s 0). The dimensionless mass and centri- 
fugal force coefficients are given by 


Mn 


M 22 


M33 



(s — 8 m ds 


( <t> m ds 


% 2 rds 



m 12 = 0 


t 12 = 0 


M13 = 


M 


23 = 



(s - 8) x m m ds 


^ f 0 m ds 


t 13 = - m 13 



(s -8) 



% x m m ds 


M ji = M ij 


T ii 


T-- 


The generalized forces are integrals over the blade span involving section lift 
and moment: 
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1 1 an ii 


1 



u = s + fJL sin i jf 

while and C m are section lift and moment coefficient, respectively. 

The computer coding used to evaluate the blade elastomechanical param- 
eters was kept separate from the main computer program. The latter program 
could, therefore, be employed without difficulty to analyze blades with differ- 
ent hub restraints than the ones treated here. 

Analysis Procedure 

At a given azimuth position tA , the section lift and moment coeffi- 
cients are computed at specified radial position, deriving section plunging 
rate, pitch angle and pitch rate from the blade response at that azimuth posi- 
tion. The integrals defining the generalized forces are then evaluated. The 
integrands of those integrals are assumed to vary linearly between the sec- 
tions at which the load is evaluated and to go to zero parabolically at the 
blade tip and the blade root. 
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The integration in azimuth angle 0 (i. e. , in time) to obtain blade 

response caused considerable difficulty in the early stages of the study. The 
complexity of the aerodynamic representation dictates that the azimuthal step 
size be as large as possible. Further, the logic involved in evaluating the 
load at stall onset and during unstall would be extremely difficult to adapt to 
an integration scheme such as the Runge-Kutta method, which requires trial 
values of the loading to be computed at intermediate azimuth angles. There- 
fore, a type of predictor -corrector scheme was devised which avoided com- 
puting the loading at a given azimuth position more than once. Specifically, 
a first estimate of the response at azimuth position 0 + A0 was obtained 
by linearly extrapolating the generalized forces, using their values at 
0 - A if/ and 0 , and integrating the equations of motion analytically 

for a linear variation of the forcing functions. The resulting response at 

+ A 0 was used to compute the section loading and corrected values for 
the generalized forces. Then, the blade response at 0 + A0 was recom- 
puted, using the corrected values for the generalized forces at 0" + A0 , 
again assuming a linear variation for the forcing functions. 

During initial calculations, the torsional degree of freedom exhibited a 
numerical instability. After trying various remedies, it was finally deter- 
mined that the last step in the integration scheme was the one causing the 
problem. When that step was simply omitted, the numerical problem was 
eliminated and results of acceptable accuracy were obtained, provided the 
integration step was no larger than about one -sixth of the period of torsional 
oscillations. This required computer running times which were considerably 
greater than had originally been anticipated, but which were still within 
practical bounds. 

For all results presented here, the integration step size was six degrees 
of azimuth, which gives a minimum of eight steps per cycle of torsional 
oscillation for the blades analyzed. Typically, one revolution of a blade 
required about 45 minutes of 360/75 computer time. 
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RESULTS OF COMPUTATIONS 


Configurations Analyzed 

For comparison with flight test data, analyses were carried out for 
both an H-34 and a CH-54B rotor blade. The basic H-34 configuration was 
used in the investigations of the effects of vibratory characteristics and of 
boundary layer control. Blade and flight parameters of the two aircraft, 
obtained primarily from Refs. 11 and 12, are listed in Table 1. Bending and 
torsion mode shapes for the H-34 are plotted in Fig. 4. The bending mode 
shape for the CH-54B is very nearly the same as that of the H-34, and the 
torsion mode shape differs only in the deflection at the root, which is . 0909, 
compared with . 0969 for the H-34. Blade bending stiffness and mass distribu- 
tions are shown in Figs. 5 and 6. 

The blades of both aircraft have an NACA 0012 section, modified by a 
small trailing -edge extension. The blade section is shown in Fig. 7, together 
with a table of offsets. The static lift and moment coefficient variations with 
angle of attack, obtained from Ref. 18, are plotted in Fig. 8, together with 
the computed variations, taken from Ref. 9. 

In the calculations, for a given advance ratio, collective pitch angle was 
first varied parametrically until the desired thrust coefficient was obtained. 
The value of Crj, used for the H-34 is . 0057 and for the CH-54B, . 010. 


Effect of Advance Ratio - Comparison with Flight Test 

H-34. — The data from Flight 15 of Ref, 11 were selected for detailed 
comparison with computed blade response and loading. During this flight, the 
blades were undergoing severe torsional oscillations due to stall. The air- 
craft was in level flight at a speed of 105 kts and an advance ratio of . 29, The 
variation with ^ of the blade torsional moment at r/R = . 15, showing 
large amplitude oscillations at between 6 and 7 cycles per revolution in the 
third and fourth quadrants, is plotted in Fig. 9. For comparison, the tor- 
sional moment for Flight 13 of Ref. 11 (level flight, V£ = 99 kts, ^ = . 27) 
is also plotted, there being no evidence of stall-induced oscillations during 
that flight. The variations of section lift and moment coefficients with ft for 
these flights, computed from the data of Ref. 11, are shown in Fig. 10 for 
r/R of . 55, . 75, . 85, and . 95. The loading variations for the higher advance 
ratio are symptomatic of dynamic stall. The lift coefficient rapidly rises to 
a value considerably in excess of the maximum static value of about 1. 2, 
while at the same time the moment coefficient drops sharply to a large nega- 
tive value. The lift coefficient then decreases while the moment coefficient 
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TABLE 1 


BLADE AND FLIGHT PARAMETERS 






BENDING STIFFNESS, N-M 



0.6 0.8 1 .0 
R 

FNESS DISTRIBUTIONS 


19 








Xj/c 

Yu/c 

0 

0 

0.0110 

0.0170 

0.0220 

0.0230 

0.0330 

0.0270 

0.0540 

0.0340 

0.0760 

0.0390 

0.1087 

0.0445 

0.1521 

0.0493 

0.2065 

0.0527 

0.2500 

0.0542 

0.3043 

0.0547 

0.3478 

0.0541 

0.4130 

0.0520 


x-|/c 

Yu/c 

0.4564 

0.0499 

0.5000 

0.0472 

0.5434 

0.0439 

0.6086 

0.0383 

0.6521 

0.0343 

0.6955 

0.0300 

0.7607 

0.0230 

0.8042 

0.0181 

0.8477 

0.0127 

0.891 1 

0.0070 

0.9346 

0.001 1 

1.000 

0.0011 


ro/c = 0.0143 


Figure 7 AIRFOIL SECTION AND OFFSETS 
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Figure 9 VARIATIONS OF BLADE TORSIONAL MOMENT AT r = 0.15 R DURING FI-34 FLIGHT TESTS 
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Figure 10 VARIATIONS OF BLADE SECTION LIFT AND MOMENT COEFFICIENTS DURING H-34 FLIGHT TESTS 













rises to a somewhat more moderate value. Unstall is marked by the return 
of the moment coefficient to near zero. 

The results of the analysis of the H-34 blade at ft = . 27 are in good 
agreement, in terms of both the extent of stall and the variation of torsional 
moment, with the flight test data. Increasing advance ratio brought about the 
expected increase in blade stall and resulting large -amplitude torsional 
oscillations, but at a somewhat higher speed than in the flight test, due to a 
more gradual growth in the stall zone with increasing ft . The computed 
blade root torsional moments for advance ratios of . 31 and . 33 as well as . 27, 
are shown in Figs. 11 and 12, together with the flight test results. There is 
some evidence of torsional oscillations in the fourth quadrant at ft = . 31 
(Fig. 11) but not as severe as what was measured. At ft = . 33 (Fig. 12), 
oscillations roughly comparable to the measured ones are seen. 

To show the relation between the severity of the oscillations and the 
extent of blade stall, the computed stall zones for advance ratios of . 31 and 
, 33 are mapped in Fig. 13, together with the measured stall zone at ft = . 29. 
The computed stall zone for the lower advance ratio is seen to be somewhat 
smaller than what was measured, with stall reaching to about r/R = . 85, 
while the computed extent of stall for ft = . 33 is nearly coincident with the 
flight test results. Thus, a relatively small increase in the radial extent of 
stall caused a large increase in the amplitude of the torsional oscillations. 

The more rapid onset of stall-induced oscillations with increasing 
seen in the flight test results could be due in part to differences in aircraft 
gross weight. There is an uncertainty of about 5% in the weight, and hence 
rotor thrust, reported in Ref. 11, so the data for Flight 15 ( ft = . 29) could 
have been taken with the rotor operating at a somewhat higher thrust level 
than for Flight 1 3 ( \i - . 27), thereby causing extensive blade stall at only 
a slightly higher advance ratio. In any case, because of the good agreement 
in the extent of stall and severity of the resulting torsional oscillations, in 
what follows the computed results for ft = . 33 will be regarded as comparable 
to the flight test measurements at ft = . 29. 

The computed variations of section lift and moment coefficients with 
if/ for ft = . 27 and ft = . 33 are shown in Fig. 14. The azimuth positions 
of stall and unstall have been marked on the abscissa. In comparing these 
results with the measured loading (Fig. 10), it is seen that the essential 
features of the unsteady stall process are reproduced. Substantial lift over- 
shoot is accompanied by a rapid drop in moment coefficient, with a partial 
increase in moment shortly after the drop in lift, followed by a return of the 
moment to near zero at unstall. 
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Figure 13 COMPARISONS OF COMPUTED AND MEASURED STALL ZONES 
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CH-54-B . — It is reported in Ref. 12 that the CH-54B encountered stall - 
induced torsional oscillations during flight tests at a speed of 115 kts, or an 
advance ratio of . 29. Analysis of the CH-54B at that advance ratio produced 
stall out to a radius of # 7 R, but there was no appreciable torsional oscillation. 
Increasing the advance ratio to . 33 caused an additional 15% of the blade to 
stall, which resulted in excessive torsional oscillations. 

The computed pushrod load variations with ft for the two advance 
ratios are shown in Fig. 17. The measured variation given in Ref. 12 is 
plotted without an ordinate scale, so a direct comparison is difficult. The 
measured oscillations are somewhat more severe than the computed ones, 
however. The former give a variation from maximum to minimum of about 
4200 lbs, while the corresponding variation from the analysis is 2300 lbs. 

Also, the measured oscillations persisted through the first quadrant of the 
rotor plane, while they are confined to the fourth quadrant in the computed 
re suit. 

The computed section lift and moment coefficients are plotted in Fig. 18 
for four radial locations. The variations with ft are seen to be quite similar 
to those obtained for the H-34. Stall progresses from root to tip, and again^ 
a small increase in the extent of stall produced a large increase in the severity 
of torsional oscillations. 


Effects of Variations in Torsion Frequency 
and Damping in Torsion 

Other investigations have provided evidence that alterations in the blade 
vibrational characteristics can alleviate the problem of stall-induced torsional 
oscillations. Dampers were installed on the rotating control linkages of a 
CH-54-B which reduced the 5th through 8th harmonics of the pushrod load by 
from 35% to 50% (Ref. 12). The analysis of Ref. 14 shows a reduction of 
about 40% in pitch-link loads would result by reducing the torsional frequency 
from 7 Q to 5. 5 0 . The H-34 rotor blade was analyzed, parametrically 
varying torsional frequency and introducing damping in torsion, to determine 
whether these effects could be detected. 

The analysis of the H-34 with p = . 33 was used as a basis for compari- 
son. Torsion frequency was changed by adjusting the control stiffness. 

Viscous damping was introduced at two tenths of critical, which is the same 
amount as was produced by the control-linkage dampers of Ref. 12. The 
damping was assumed to act without altering the undamped torsional mode 
shape. The blade root torsional moment for each case was harmonically 
analyzed, obtaining the zeroth through ninth harmonics of the rotational fre- 
quency to provide a quantitative basis for comparison. 





The effect of varying torsional natural frequency is shown in Fig. 19, 
which compares the amplitudes of the zero^ through ninth harmonics of the 
torsional moment for torsional natural frequencies (rotating) of 5, 6 0 , 

6. 6 0 , and 7. 6 Q . Little effect is seen through the fourth harmonic. 

Increasing the frequency to 7.6 fi lowered the fifth harmonic by 30%, but 
raised the sixth about the same amount and left the seventh through ninth 
nearly unchanged. Lowering the frequency to 5. 6 £2 caused about a 30% 
increase in the fifth harmonic, but significantly lowered the sixth, seventh, 
and eighth. The largest of these, the seventh, was reduced by 60%. Lowering 
the torsional natural frequency does have a significant favorable effect, then, 
in agreement with what was found in Ref. 14. 

The results of the harmonic analyses of the torsional moment for the 
three frequencies with and without damping in torsion are shown in Figs. 20, 
21, and 22. Damping caused significant reductions for all three frequencies, 
and particularly so for a frequency of 6. 6 12 (Fig. 21), At that frequency, 
the seventh through ninth harmonics are reduced to almost negligible levels, 
and the zero^, first and second are substantially reduced as well. The first 
harmonic, which is the largest, was reduced by 26%. These results indicate, 
then, that the introduction of damping in torsion can significantly reduce 
stall-induced torsional oscillations, which is in concurrence with the findings 
of Ref. 12. 


Investigation of Boundary Layer Control 

A preliminary study of the effectiveness and feasibility of using 
boundary layer control to reduce stall -induced torsional oscillations was 
conducted. To determine effectiveness, stall was not allowed to occur over 
a specified blade segment during the analysis of blade loading and response. 
Estimates were then made of pumping rates and power requirement s, using 
available experimental data as a basis. 

The H-34 rotor blade at an advance ratio of . 33 was again analyzed. 
The procedure for evaluating the spanwise loading distribution was altered 
to account for boundary layer control, as follows. If the segment of blade 
designated to be controlled extends from r = to r = rg, then three of the 
six radial points chosen for evaluating section loading were r a , rg, and the 
midpoint of the controlled segment, r c = (r^\ 4- rg)/2. It was assumed that 
the loading on the controlled segment was constant and equal to its value at 
r = r c . Spanwise integrations were then divided into three parts, one from 
the blade root to r^, another from r^ to rg, and the third from r-g to 
the tip, in order to allow for the possible loading discontinuities at r^ and/ 
or t - q . Stall was not allowed to occur at r = r c , but the leading-edge 
bubble was still analyzed there. The results of those analyses were used in 
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Figure 18 COMPUTED VARIATIONS OF CH-54B BLADE SECTION 
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the estimates of pumping and power requirements. The blade root torsional 
moments were again harmonically analyzed for comparison with the case of 
no boundary layer control. 

The blade was analyzed for lengths of the controlled segment of blade, 
l c , of . 05 R, .1 R, and . 15 R. The results are presented in Figs, 23, 24, 
and 25, where the amplitudes of the fifth through eighth harmonics of the tor- 
sional moment are compared for different values of r c# The lower harmonics 
were generally reduced by boundary layer control, but the reductions are not 
substantial and no clear trends are evident, so only the higher harmonics are 
shown. 


As one would expect, increasing l c reduces the severity of the tor- 
sional oscillations. For l c of . 1 R and . 15 R (Figs. 24 and 25), it is appar- 
ent, as well, that control near the tip is more effective than control over an 
inboard section. The greatest overall reduction in the amplitudes of the 
higher harmonics was obtained with l c = , 15 R and r c = . 9 R, with the largest 
of these, the seventh, reduced by 95%. The variation of the torsional moment 
with \jj for that case is compared with the result for no control in Fig. 26. 


In order to estimate pumping and power requirements, data were used 
from the wind tunnel tests on an airfoil with 9% thickness ratio at a Reynolds 
number of 1. 15 x 10^ reported in Ref. 20. Boundary layer control in those 
tests was effected by suction over a surface area near the leading edge com- 
prising 2. 75% of chord. The results are summarized in Fig. 27, taken from 
Ref. 20, which shows the lift coefficient vs. incidence for various suction 
rates. The suction rate required to just prevent stall is see n to increase 
approximately linearly with incidence, with a value for Cq \/ Re of 3. 61 
being sufficient to increase the stall angle from 10 degrees to 16 degrees; the 
suction coefficient Cq = Q/(Uc), Q being th e vo lume flow per unit span and 
U the free-stream speed. The product Cq \/Re is used because results are 
nearly independent of Re for a given value of that parameter. 


In order to compute Q from these results, the equivalent of airfoil 
static angle of attack is needed. As is apparent from the results presented 
previously, the actual angle of incidence is not a good measure of how far the 
airfoil is above the stall condition in unsteady flow, so. results from the 
analysis of the leading-edge bubble were used. Bet A f denote the required 
increase in pressure coefficient across the bubble, and the increase 

which is possible, as defined in Ref. 8. Then, a parameter £ c is defined by 


- 



\ > A p 
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Figure 24 HIGHER HARMONICS OF TORSIONAL MOMENT WITH £ C /R = 0.1 







Figure 27 VARIATION OF LIFT COEFFICI ENT WITH INCIDENCE AND 
SUCTION QUANTITY (FIG. 9 OF REF. 20) 




to provide a measure of how far above stall an airfoil is under either steady 
or unsteady conditions. To relate £ c to static angle of attack, the leading- 
edge bubble was analyzed in steady flow with stall precluded for several 
angles of attack above the stall incidence, with the result shown in Fig. 28. 
The parameter f c was defined so it would vary linearly with angle of 
attack near the static stall angle, as can be seen is the case from Fig. 28. 

The data of Fig. 27 were combined with the results of Fig. 28 to give 

the plot of Cq \/ Re vs. shown in Fig. 29. The maximum value of £ c 

at r c for each run was determined and the corresponding value of Cq /Re , 
taken from Fig. 29, was used to compute required suction rate. 

The power required was computed as the product of the total volume 
flow per blade, QtOT> the increase in total pressure supplied by the pump, 
and the number of blades, Ng. Neglecting duct losses and assuming exhaust 
is to atmospheric pressure, the total required power is given by 


p s - — n b Qtc>t p (% 2 + ^ 2 r c 2 ) 

where q g is the magnitude of fluid velocity at the edge of the boundary layer 
at separation. 

Results are summarized in Table 2. Calculations for r c /R = . 95 and 
1 C /R = . 05 could not be performed because A r never exceeded A p during 
the run but, as might be expected, the reduction in torsional response was 
not appreciable for that case (see Fig. 23). The suction rates and power re- 
quired for all cases are relatively modest. It should be noted that require- 
ments do not necessarily increase with effectiveness, since the case resulting 
in the most reduction in torsional oscillations (r/R = . 9, 1 C /R = . 15) does 
not require the most power of the cases analyzed, by a considerable margin. 

In light of these results, the primary factor in determining whether boundary 
layer control is implemented would appear to be the mechanical complexity 
and resulting cost, rather than duct size or power requirements. 
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TABLE 2 

SUCTION RATE AND POWER REQUIREMENTS 


r c /R 

V R 

c max 

Cq \/Re 

10" 6 Re 

q s /( fl R) 

Qtot m ^/ sec 

P s hp 

. 75 

. 05 

. 527 

4. 95 

1.47 

1.42 

. 0578 

17. 7 

. 85 

. 05 

. 405 

2. 18 

1. 84 

1. 65 

. 0286 

11. 7 

. 95 

. 05 

-- 


-- 

-- 


-- 

. 75 

. 10 

. 553 

5. 90 

1.49 

1. 50 

. 139 

46. 5 

. 85 

. 10 

. 369 

1. 97 

1. 82 

1. 69 

. 0513 

21. 7 

. 95 

. 10 

. 233 

1. 25 

2. 23 

1. 99 

. 0360 

20. 7 

. 7 

. 15 

. 596 

7. 50 

1. 35 

1. 39 

. 252 

72. 3 

. 9 

• 15 

. 283 

1. 50 

2. 00 

1. 83 

.0614 

30. 3 






CONCLUSIONS 


An analysis of helicopter rotor blade torsional oscillations due to stall 
has been carried out, using an analytic dynamic stall representation to 
determine aerodynamic loading. The following conclusions can be drawn from 
the results obtained. 

1. The analytic representations employed will predict the 
torsional oscillations due to stall and their relationship 
to blade and flight parameters. In a comparison of 
data from flight tests with results of analyses, the 
severity of the oscillations and their relation to aircraft 
forward speed agreed fairly well, while the radial and 
azimuthal extents of stall and the azimuthal variations 
of section loading were in good agreement. 

2. The amplitudes of the higher harmonics of torsional 
oscillations can be significantly reduced by either re- 
ducing the torsional natural frequency or introducing 
viscous damping in the torsional degree of freedom. 

3. Boundary layer control would be an effective means for 
reducing the higher harmonics of torsional oscillations 
due to stall. Its implementation would not require 
excessive power or suction rates. 
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